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Abstract. 

Gravitational lenses on galaxy scales are plausibly modelled as having ellipsoidal 
symmetry and a universal dark matter density profile, with a Sersic profile to describe 
the distribution of baryonic matter. Predicting all lensing effects requires knowledge 
of the total lens potential: in this work we give analytic forms for that of the above 
hybrid model. Emphasising that complex lens potentials can be constructed from 
simpler components in linear combination, we provide a recipe for attaining elliptical 
symmetry in either projected mass or lens potential. We also provide analytic formulae 
for the lens potentials of Sersic profiles for integer and half-integer index. We then 
present formulae describing the gravitational lensing effects due to smoothly-truncated 
universal density profiles in cold dark matter model. For our isolated haloes the density 
profile falls off as radius to the minus fifth or seventh power beyond the tidal radius, 
functional forms that allow all orders of lens potential derivatives to be calculated 
analytically, while ensuring a non-divergent total mass. We show how the observables 
predicted by this profile differ from that of the original infinite-mass NFW profile. 
Expressions for the gravitational flexion are highlighted. We show how decreasing 
the tidal radius allows stripped haloes to be modelled, providing a framework for a 
fuller investigation of dark matter substructure in galaxies and clusters. Finally we 
remark on the need for finite mass halo profiles when doing cosmological ray-tracing 
simulations, and the need for readily-calculable higher order derivatives of the lens 
potential when studying catastrophes in strong lenses. 
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1. Introduction 

A large number of cold dark matter (CDM) N-body simulations agree that the haloes 
formed have, on average, a universal broken power law density profile. While there is 
some debate over the logarithmic slope of the profile within the break, or scale radius r s , 
and there is some scatter between haloes, it seems that the original profile of [l| (NFW) 
still provides a reasonable fit to the simulation data. This profile seems to be a generic 
feature of haloes formed in the hierarchical model of structure formation. 

A number of authors have made significant progress in understanding non-linear 
effects in structure formation using halo models, where the number density, correlation 
function and mass density profile of CDM haloes are fitted to the numerical simulations 
and then used in a simplified model of large scale structure 0, 0, 0]. The NFW profile 
has played a prominent role in this enterprise. One application of the halo model is in 
the investigation of the halo occupation distribution, that characterises the substructure 
within a larger halo. This has long been one of the more controversial topics in CDM 
theory, with predictions and observations often at odds. In order to build up an accurate 
picture of a hierarchical mass distribution, the stripping of the sub-haloes by tidal 
gravitational forces must be modelled j^l, @]; measurements of halo stripping form an 
important test of the detailed predictions of the CDM simulations 0, 0] • 

Moreover, it is now clear that the effect of baryons on the shapes and profiles of total 
mass distributions cannot be ignored. In galaxies, the stellar component of the mass 
dominates at small radii giving rise to a peakier observed total density profile than 
seen in pure dark matter simulations [9]. The surface brightness profiles of massive 
galaxies seem to be consistently well- fitted by a Sersic profile of index 3 — 4 [10 1; a 
logarithmic slope of —2 in total density in the inner regions appears ubiquitous 111 ]. 
Moreover, the dark matter profile itself is expected to steepen during the formation of 



the galaxy, by the process of adiabatic contraction 12J, ll3j]. Typically this leads to more 
centrally concentrated, rounder haloes [hi, [l5j]. The details of the mass distributions 
of real galaxies are therefore a probe of the galaxy formation physics claimed by the 
simulations. 

Gravitational lensing allows us to probe the mass distributions of galaxies, groups 
and clusters in a unique way. Insensitive to the dynamical state of the lens system, both 
weak and strong lensing effects depend only on the projected (and scaled) gravitational 
potential. Gravitational lensing has already been used to investigate the density profile 



in galaxy clusters [ly, |l7|, |l8|, |l9j. Substructure studies have also been undertaken, 



making use of the galaxy-galaxy lensing effect in clusters 20|. T he g alaxy scale halo 
mass profiles have also been measured, using both strong lensing 2ll . |22j, |9J] and, in a 
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more statistical fashion, galaxy-galaxy weak lensing 

Lensing studies provide direct tests of the CDM simulations, and typically involve 
(at some point) fitting the parameters of an NFW-like model to the data. However, 
this model is also well-suited to a more general analysis, building up a data model from 
a linear combination of NFW-like potentials (25)]. This approach has applications in 
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substructure characterisation, and also template-based cluster finding. Characterising 
stripped substructure both require an accurate treatment of the outer regions of haloes. 
However, in order to measure accurately density profile slopes and concentrations, the 
baryonic mass component must be included. 

In the perpetually applicable thin lens approximation it is the projected Newtonian 
gravitational potential that gives rise to the gravitational lensing observables. Making 
the simplifying assumption that projected stellar mass density is proportional to optical 
surface brightness leads us to seek the potential that corresponds to the Sersic density 
profile. Likewise, for the dark matter component we require a lens potential that 
corresponds to the universal profile seen in simulations, but that also includes the effects 
of tidal stripping. 

Analytic lens potentials are convenient to work with: they, and their derivatives 
that are needed for lensing data modelling, can be computed quickly and accurately; the 
introduction of ellipticity to the halo can be done very straightforwardly; more complex 
potentials can be constructed by simple linear combination of analytic functions. This 
last feature allows concave isodensity contours to be avoided in the case of high ellipticity. 
It also allows total density distributions to be constructed from mixtures of dark and 
luminous matter. 

In this work we present analytic forms for a smoothly truncated universal CDM 
gravitational potential, and also for the potentials corresponding to a subset of the 
Sersic profiles. If the underlying potential is analytic, so are all the derivatives needed 
in gravitational lens studies. An outline of the paper is as follows. In Sections [2] and [3] 
we present our suggested analytic potential models, for both dark and baryonic matter 
components, and outline the derivation of the quantities relevant to gravitational lensing. 
We leave the full formulae to an appendix, but in Section @] we plot the predicted 
observables, and compare them to those from an unstripped baryon-free NFW form. In 
Section [5] we briefly discuss our results, and point the reader towards some publically 
available computer code. 



2. Smoothly truncated dark matter haloes 

The NFW profile for the CDM density p of a halo is 
/ \ _ 4(5 c p crit _ M 



The characteristic overdensity 5 C is the density at the scale radius r s , in units of the 
critical density p cr i t . Alternatively, we can express this as p(r s ) = 5 c p C r\t = M /(167ir^). 
The NFW profile is analytically integrable along the line of sight; the most frequently 



used formulae for the weak lensing shear 26( and strong lensing image positions 27 



were 



derived assuming the integral to extend over all space. Given that the NFW profile has 
divergent total mass, 0] suggest using a modified form that is sharply truncated at the 
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virial radius. The projection integral is then more realistic, with only mass physically 
associated with a finite-sized halo being modelled. 

We might expect real CDM haloes to be truncated due to tidal effects; a step- 
function density cutoff may not offer a very physical picture of the edges of haloes. 
With the 0] mass distribution, the lensing deflection angle and shear are tractable 
(if somewhat less simple), and the actual potential is worse still, involving (at least) 
polylogarithms. Also, the convergence and shear are not differentiable at the truncation 
radius. A power-law cutoff in the potential is more attractive in this regard. We should 
insist that the truncated profile match that of NFW as closely as possible within the tidal 
radius, which is introduced as the third parameter of the profile. This is important for 
the results from previous work on fitting the outputs from N-body simulations pertaining 
not only to the density profile but also the mass function [281 ]. 

With these desiderata in mind, we suggest the following functional form for a 
smoothly truncated universal 3-d mass density profile: 

p(r) = ^ = M ° (^L\ n (2) 

PU '^)(l + ^) 2 (l+(^)y 4 ^(r + r s ) 2 ^ 2 + nV ' U 
Here, r t is a new parameter which should correspond to the tidal radius for tidally 



truncated halos [29|. The parameter n controls the sharpness of the truncation; we 
will investigate the cases n = 1,2. For relatively isolated haloes, we expect the tidal 
radius to be much larger than the scale radius. We define r = r t /r s , expecting r >> 1. 
Note that r is not necessarily the "concentration" parameter, defined as the ratio of a 
"virial" radius to r s . In the left-hand panel of Figured! we plot this profile in the usual 
way (with logarithmic axes), and compare with the original NFW profile. We show 
the effect of the tidal radius in providing a smooth edge to the halo. In the right-hand 
panel of Figure [1] we show the integrated projected mass profile, for the same set of 
density profiles. In this panel radius is projected radius, and the logarithmic divergence 
of the original NFW profile can be clearly seen. Projected mass within some appropriate 
radius is (approximately) the quantity that is best constrained by gravitational lensing 
- in Section H] we show the predicted observables of gravitational lensing in more detail. 

After some experimentation we found that the form recommended above is indeed 
the simplest one that gives an analytic potential while ensuring a non-diverging total 
mass for all values of the tidal radius. In the case where the tidal radius rt is outside the 
scale radius r s , the n = 1 profile falls off as r~ 5 , steep enough to mimic a sharp cutoff. 
If the tidal radius were to lie inside the scale radius, then the density would decrease as 
r~ 3 initially before turning over to r~ 5 outside r s . Since this would imply some memory 
of the original halo after the presumably violent act of tidal stripping, we suggest that 
if r t < r s , the n = 2 version of the density profile be used. For r t < r s , this profile 
turns over to r~ 5 at r t , which is effectively a sharp cutoff. The further turnover to r~ 7 
at r s > r t has little effect. 

The close agreement of the unstripped halo with the original NFW profile is 
comforting. For example, for a halo with a concentration of 10, if we set r t to twice 
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Figure 1. Density (left) and integrated projected mass (right) profiles for NFW haloes 
with various truncation schemes. The solid lines indicate the original NFW halo, with 
and without a hard cutoff at r = rt/r s — 10. Dotted lines show the n — 1 cutoff 
prescription, with r — 10, 20. Dashed lines show the n = 2 cutoff prescription, again 
with r = f 0, 20. For smoothly truncated models, ratios of masses outside the virial 
radius to the virial mass, [M tot - M (< 10)]/M(< 10), are 17% for (n, r)=(l, 10), 4.6% 
for (2, fO), 36% for (1, 20), and 17% for (2,20). Note that this ratio is infinity for the 
original NFW halo because it has divergent total mass. 



the virial radius (r = 20), the masses contained within the virial radius are the same 
to within 6%. Of course the total mass of the truncated halo is 50% larger than the 
virial mass, while the total mass of the untruncated halo (formally) diverges. We will 
thus take r t = 2r 2 oo as our fiducial tidal radius for an unstripped halo. We note that 
this choice of the truncation radius is simply a working assumption in this paper, and 
the more appropriate value should eventually be obtained in iV-body simulations and 
observations. 



3. Lensing by stellar mass in galaxies 



The Sersic profile, found to fit well the optical surface brightness I of undisturbed 
galaxies 



30 



is 



Kr) 



exp 



r 



l/n 



(3) 



where the effective radius r P is the radius within which half the flux is contained, and n 



is the Sersic index. For elliptical galaxies, an index of around 4 is often seen 31] . while 
the characteristic exponential profile of galaxy disks corresponds to a Sersic index of 1. 
In fact, a broad range of Sersic index values have been seen in fits to observed galaxy 



light profiles 32j. In addition, there has been some arguments that density profiles of 



CDM haloes can also be fitted well by the Sersic profile with an index of 2 — 3. [33j, [34 
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Assuming that stellar mass follows light, we can substitute surface mass density E 



for surface brightness in equation [3 



In the appendix, we show that the lens potential 



sourced by this mass distribution [35j is analytically tractable for integer and half-integer 



n. 



4. Predicted observables 



As we show in the appendix, both density profiles introduced above (truncated NFW 
and Sersic) have analytic lens potentials. (The NFW profile has an analytic three- 
dimensional potential, which can itself be projected analytically.) The expressions for 
the lensing potential, while somewhat lengthy, are rapidly calculated and differentiable 
to all orders. In this section we plot some of these derivatives, pointing out their 
application in gravitational lens data modelling. We note that in the limit of radii 
beyond the tidal radius the lensing properties of our model haloes do indeed approach 
those of a point mass, as required. 



4-1. Weak lensing 

We first address the issue of not truncating the NFW profile on the weak lensing shear 



and convergence (see [36|] for a good introduction to these quantities). The lefthand 
panel of Figure [2] shows the convergence (projected mass density) profile for the set of 
haloes first introduced in Figure [TJ 

We see that using a truncated profile with the same virial mass somewhat reduces 
the predicted lensing effect. The corresponding shear profiles are plotted in the right- 
hand panel of Figure [3 Taking the central density profiles to be equal (mimicking a 
well-constrained central strong lensing region, for example), for r = 20 we note that the 
virial mass (mass within 10r s ) is 6% less for the n = 1 truncated halo. The shear for the 
two profiles only differs by 3%, however. The total projected mass within the projected 
virial radius is some 12% lower than that of the untruncated profile. Lastly, the surface 
density of the truncated halo is 30% lower at 10r s . 

The difference in reduced shear — k\ thus depends on the absolute value of 
the convergence, relative to the critical surface density, but can be significant. Very 
roughly, from Figure [2] we expect that different truncations examined in this paper can 
yield ~ 10% difference in 7 at around the virial radius. Although this is smaller than 



the accuracy of shear measurements for most massive clusters of galaxies (e.g., [37|), the 



accuracy can be reachable in the weak lensing analysis of stacked cluster samples (e.g., 



4-2. Strong and intermediate lensing 

Figure [3] shows the amplitude of the deflection angle for a strongly lensed source. Here 
we see that using a truncated profile has very little effect on the deflection angle in 
the regime where it is measurable as a multiple-image separation (r < r s ). This is 
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radius x radius x 



Figure 2. Convergence (left) and shear (right) for NFW haloes. The curve line 
styles are the same as in Figure [TJ Projected mass density £ (directly proportional to 
convergence k) is plotted on the left, showing that a hard cutoff in density results in a 
softer, but still non-differentiable, cutoff in the convergence. Actually plotted on the 
right is (£) — S, which is directly proportional to the shear 7 for axisymmetric haloes. 
Notice that a hard cutoff in density means the shear is finite, but not differentiable at 
the cutoff radius. 
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Figure 3. Deflection angle (left) and flexion (right) for NFW haloes. Curves are the 
same as in Figure [TJ Notice that the hard cutoff in density causes the flexion to diverge 
at the cutoff radius. As the flexion approaches —00, so it changes sign as well. This 
can be seen in Figure [21 where the shear actually starts to increase as the cutoff radius 
is approached. 



unsurprising given that the strong lensing is dominated by the central part of the profile 
which is, by design, little changed in our new model. 

In contrast, the so-called "flexion" may be more strongly affected by the truncation 
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because it is essentially the higher-order derivative of the lens potential. In addition it 
is measurable over a wide range of scales from the Einstein radius to (in a statistical 
sense) the virial radius. In Figure [3] we plot the third derivative of the lens potential: 
the most interesting component of flexion for circular symmetry is in fact just this, the 
radial gradient of the shear. Marked differences in signal strength arise when the haloes 
are truncated. 

A plausible model for an elliptical galaxy lens consists of two parts: the stellar 
component and the dark halo. Modelling the stellar component as an n = 4 Sersic 
profile and the dark halo as a truncated NFW profile with concentration 10 and r = 20, 
a reasonable fit to lensing data can be made |39j. The salient feature is that the total 
mass profile is approximately isothermal. This can be arranged with the following 
prescription (for r = 20 and n = 4): 

^™«4.75^. (4) 
M deV r e 

The NFW mass is the total mass. The virial mass (mass within 10 scale radii in this 
case) is 0.66 of the total. The broad isothermal region obtained by this prescription 
is illustrated in Figure HI We note that adiabatic contraction 12j, ll3j, ll4j modifies the 
NFW profile, leading to more centrally concentrated profile of dark matter. However, 
the effect of adiabatic contraction is most pronounced at the very center of the halo 
where the baryonic (Sersic) component is dominated, and thus its effect on the total 
mass profile is not substantial. 

Finally we move to two dimensions and illustrate the construction of elliptically 
symmetric isophotes in the convergence distribution. It has been noted 40j that the 
isodensity contours of an elliptical lens potential can become dumbbell-shaped at low 
values of the axis ratio. In the appendix we show how isodensity contours that are 
elliptical to third order in the axis ratio can be constructed following a simple recipe. 
In Figure [5] we illustrate this procedure, showing the constituent concentric elliptical 
potentials and the resulting convergence contours. It is found that our new model not 
only avoids the unphysical concave isodensity but also gives much better fit to the 
elliptical isodensity than the elliptical lens potential. For large axis ratios (3:1 is as 
large as is feasible by our technique), the isodensity contours are slightly disky. We note 
that this recipe preserves the radial profile of the (self-similar) constituent potentials. 
Although the procedure is derived assuming a pure power-law, the nearly isothermal 
distribution of the composite model illustrated in Figure H] suggests that our prescription 
is useful for such composite model, at least as long as the ellipticities of Sersic and NFW 
components are similar. 



5. Discussion 



We have introduced a simple smoothly-truncated extension of the NFW density profile. 
To date the majority of cluster and galaxy lens modelling that has been performed using 
the NFW profile has used the untruncated profile. We find that, if haloes are indeed 



NFW lenses 



9 



10° 



10- 



10- 



10- 



10° 10 1 
x (half light radii) 



10 2 



10- 



10° 10 1 
x (half light radii) 



10 2 



Figure 4. Sersic (n=4) profile combined with NFW profile. We plot the convergence 
(left) and circular velocity squared M(x)/x (right) associated with each component 
(there are two curves for the NFW profile: one untruncated and one with r = 20), 
along with the total (dashed line). The truncated NFW profile is 70 times more massive 
than the Sersic profile (the virial mass is 50 times larger than the stellar mass), and 
the scale radius is 15 times larger than the half light radius. With these reasonable 
parameters, it is clear that the total profile is nearly isothermal (logarithmic slope of 
-1 in convergence) around the half light radius. In fact, we find that (for r = 20), the 
relation M^pw/M^cV = 4.75 r s /r e gives a flat region in velocity dispersion. 




Figure 5. Convergence (isodensity) contours for three stacked elliptical potentials. 
True ellipses are shown as solid curves, three stacked potentials are shown as dashed 
curves, and single elliptical potentials are shown as dotted curves. The left panel 
illustrates the case ej so = 0.6 (axis ratio 2:1), while the right panel illustrates e; so = 0.8 
(axis ratio 3:1). The slight diskiness of the isodensity contours can be seen especially 
for the eiso = 0.8 case. The fitting procedure is described in Appendix C. 



NFW lenses 



10 



tidally-moulded leading to the kinds of smooth truncation that we propose, then the 
masses of the haloes may have been overestimated by some 10% or so during a weak 
lensing analysis. This number pertains to the situation where the inner profile is inferred 
to be the same for each model profile, as might be the case when good strong lensing 
data are available. We find that the smooth truncation of a halo does not significantly 
affect the deflection angles at the image positions (which lie typically well within the 
scale radius). If strong lensing data are not available then the degeneracy between 
the truncation radius and the halo mass will give rise to a broader inferred marginal 
probability distribution for the halo mass, with the mean shifting to lower values than 
for the untruncated profile. 

The truncation of galaxy haloes in field galaxy-galaxy lensing is always likely to 
be masked by the effects of large scale structure on the outer parts of the mass profile 
(the "two-halo" term). However, in clusters a 10% systematic error is comparable to 
that introduced in other parts of a current weak (plus strong) lensing analysis. The 
uncertain background galaxy redshift distribution, additional mass along the line of 
sight, cluster member galaxy contamination, projection effects, and shear calibration 
errors can easily be of order 10% in the halo mass. However, as survey sizes increase, 
and the goals of cosmological cluster-counting experiments become loftier, uncertainties 
such as that introduced by halo truncation may become important. 

Sharp truncation [4] introduces discontinuities in the shear that are unlikely to 
cause problems in data modelling; however, the same may not be true about flexion, 
where a singularity appears at the truncation radius. Smooth truncation (or indeed, no 
truncation) avoids this problem. 

Galaxy-galaxy weak lensing studies within clusters of g alaxies have already 

The combined Sersic 



41 



succeeded in producing a measurement of a truncation radius 
plus NFW model currently popular in field galaxy lens modelling could profitably be 
applied in a cluster galaxy-galaxy lensing. The photometry provides extra constraints 
on the stellar mass part of the density profile 42|, allowing the dark matter structure 
of galaxies in clusters to be probed. The model suggested here would straightforwardly 
allow strong and intermediate lensing effects to be incorporated; we are not far from 
possessing a useful sample of strong gravitational lenses lying in clusters. 

The question of how best to model tidal-truncation of dark matter haloes has been 
approached here in a phenomenological and pragmatic way: we wanted an analytic form 
for the lens potential. We believe the forms presented here would provide good fits to the 
haloes seen in numerical simulations, based on the successes of others with very similar 
profiles 0, H]. We have shown that plausible models of gravitational lenses can be 
constructed from the superposition of simple analytic profiles, including the generation 
of elliptically-symmetric isodensity contours. An interesting extension of this would be 
to attempt to build up still more complex mass distributions, from misaligned and offset 



building blocks [25|. Again, whether the haloes and sub-haloes observed in numerical 
simulations can be well-enough approximated by such a model remains to be seen. At 
present it seems that the signal-to-noise in Einstein rings is sufficiently high to constrain 
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such a more complex model |43j, [9J. 

Finally, we discuss the importance of truncation, and a smooth one at that, in when 
simulating lensing effects in large surveys. Gravitational lensing, weak, intermediate 
and strong, may be expected to be an important component of multiple-pronged dark 
energy investigation. Simulations of large fields will play a vital role in improving 
our understanding of the astrophysical systematic errors present in cosmic shear, and 
cluster mass function, measurements. Halo models are a cheap and efficient way of 
doing this, capturing the pertinent physical effects without the need for further CPU- 
expensive N-body simulations. However, ray-tracing through halo models does present 
some technical challenges 0, 44]. The smooth analytic truncation proposed here allows 
the mass budget to be balanced, while allowing all gravitational lensing effects to be 
calculated rapidly to machine precision. In fact, it has been shown that the shear angular 
correlation function becomes ~ 20% smaller if the truncation at around the virial radius 
is included, and that the calculation with the truncation shows better agreement with 
iV-body simulations [4|. 

A by-product of the future large optical and radio imaging surveys will be an 
interesting sample of strong lenses showing higher-order catastrophes beyond the usual 
cusps and folds. These systems provide very high magnifications, and are very sensitive 
to the mass structure in the lens, and as such promise to be interesting laboratories. 
However, modelling them will require an accurate multi-scale approach; we leave the 
development of this project to further work. 

The code used in this work is plain ISO C99 and can be freely downloaded from 



http : / /kipac . Stanford . edu/ collab/ research/lensing/ample/ 
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Appendix A. Truncated NFW Profile 



The Navarro, Frenk and White (NFW) profile is given by 

Mn 1 
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Defining x — r/r s , 

p (x) = - ° - ■ , 2 . A.2 

This profile describes dark matter haloes well, out to the virial radius. However it suffers 
from the deficiency that it has infinite mass. The truncation radius is defined to be a 
factor of r larger than the scale radius. With the sole motivation of allowing "simple" 
analytic forms, we propose the following truncated NF W profile: 

( \ M ° 1 ^ /A 

PT{X) ~ 47rr3x(l+x)2r2 + x2- ^ 
This form is quite similar to the NF W profile for x < r, i.e. inside the truncation radius. 
Furthermore, the total mass is finite, 

M = M (r2 ^ 1)2 [(r 2 ~ 1) Inr + ttt - (r 2 + 1)] . (A.4) 

In the limit r — > oo, we recover the logarithmically divergent mass of the NFW profile. 

For purposes of gravitational lensing, we are interested in the projected mass 
density. We first define the function 

i-./ \ cos _1 (l/x) , . . 

F(x) = r^- - (A.5) 
— 1 

This function is well defined everywhere: taking the appropriate limit F(l) = 1, and 
for x < 1, where both the numerator and denominator are purely imaginary, we choose 
the branch where F(x) > 0. Note that ArcCos from Mathematica and cacos from C99 
disagree on the sign of cos^ 1 (x) when x > 1. Note that in the x — > limit it reduces to 
F(x) = ln(2/x). We also define the following logarithm, which will appear many times, 

^ - b hw^wr) ■ (A - 6 » 

With these definitions, the projected surface mass density is given by 

S(x) =r a f d£ p T (Vi 2 + xA (A.7) 

J —CO 

M ° 71 (r2 + 1 -[l-F(x)] + 2F(x) 



r 2 s 2tt(t 2 + l) 2 [x 2 

K + L(.v)\. (A.8) 



\/t 2 + x 2 t\/t 2 + X 2 
Notice that for x — * 1, the first term requires that a limit be taken. As before, in the 
t — > oo limit we recover the result for the NFW profile, 

M 1 - F(x) „ /lnr\ . , . 

We can derive the convergence k simply by taking k = S/S crit . We will also need the 
total projected mass inside radius x, 

M proj (x) = r 2 s ! dx' 2irx"E(x') (A. 10) 

Jo 
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= Mo V + 1 + 2(* 2 - 1)] F(x) + 

+ (r 2 - l)lnr + Vr 2 + x 2 



TIT 



-71 + 



r 2 -l 



We again recover the NFW result in the r — > oo limit, 

In r 



M proj ( a :) = Mo(F(x)+ln|)+0 



(A.ll) 



(A.12) 



A crucial quantity is the shear of the gravitational field. The simplest way to derive it 
for a circularly symmetric lens is to use the mean projected surface density inside radius 
x, which is simply S = M pTO j/nr 2 . Defining r = £ — £, the shear is then 7 = r/S crit . 

We can now derive the lensing potential. Including the famous factor of two, and 
defining u = x 2 , 

s AG dx' , 2G r du' r , r-, , A , , 

= ~T ~T M P r oi (x') = — — M proj (v^). (A.13) 
c j x c j a 

The potential is thus 

, , . 2GM 1 f 2 T ro , ( ro M 

v( u > = — ^2 — ( 2 + l) 2 \ 71 T ~ ^ + u + t in [t + Vr + u J 

+ 2(r 2 - l)rv / ^ 2 T^L( v ^) + r 2 (r 2 - 1) L 2 {^i) 
+ Ar 2 {u - 1)F(V«) + r 2 (r 2 - 1) ^cos" 1 
+ r 2 [(r 2 -l)lnr-r 2 -l] \nu 

-t 2 [(r 2 - l)lnr ln(4r) + 2 ln(r/2) - 2r(r - vr) ln(2r)] |(A.14) 

The resulting potential has the correct behavior in the r — > 00 limit for all values of w, 

2 



2GM ( 







cos 1 ) br 



+ 



In 



(A.15) 



We also verify that for u » r, the potential is that of a point mass of the correct total 
mass (up to an irrelevant constant A), 

i>{u) =A + ^ [(r 2 - 1) lnr + rvr - (r 2 + 1)] In yfr 

2 s 



+ o( T - 

AGM 



= A + 



In y/u + O 



We can now calculate the derivatives of ip. The first is obvious, 
2G M proj (v^) 



u 



(A.16) 
(A.17) 

(A.18) 
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ip"(u) 



2G 



M proj (v^) 7rr s 2 S(v^) 



2GM 



v. 



u 



2nGr 2 T 
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(A.19) 



c 2 2(r 2 + l) 2 u 2 _ 
(t 2 -1)(m + 2t 2 ) 
t\Jt 2 + u 
(r 2 + 1)« 



2(1 - u - r 2 ) F(y/u) - 2{r 2 - 1) lnr 
7t(v / t 2_ +m — r) 2 



u — 



f [l - 



Again, the NFW result appears in the infinite r limit 

2GM 



u + (2 - 3w)F(v^) _ J_ 
2u 2 {u - 1) u 2 n 



+ 0(i 



(A.20) 



(A.21) 



At this point, we can calculate the time delay, deflection angle, convergence, shear 
and magnification of the truncated NFW lens. We will calculate one more derivative 
deriving the so-called "flexion" (e.g., [45l. |46I|). 

2G 



ip (u) 



2M proj (y^) _ 2 7 rr 2 S( v ^) 7rr 2 E'(V^) 



u:' 



2m 3 / 2 



2ttGV 2 
c 2 u 2 



(uV - r) 



(A.22) 
(A.23) 



2GM n 



3(r 2 + Ifu 



c 2 2(t 2 + 1) 2 u 3 {2(t 2 + u)(u-1) 2 



[uF(y/u) - 1] 



+ 



u 



2(t 2 + u)(u- 1) 
-(t 2 + 1)((m-1) + 2t 2 + 3 
2(r 2 + 1) 



[(r 2 + 1 - 2(u - l))uF(Vu) 



+ u _ 1 [F(V^) -u]+2 [2{u - 1) + 3r 2 + l] F(y/u) 



+ 



r 2 -l 



2r(r 2 + ufl 2 



[3(3r 2 + n)(r 2 + n)-r 4 ] L(y/u) 



+ 4(r 2 - 1) In 



T 



2tt(Vt 2 + u- rf 



iru 



\Jt 2 + u 

Expanding about r = oo, we find the NFW result, 



2 ( r 2 +M )3/2 



(A.24) 



(15m 2 - 20m + %)F(y/u) + 3u - 6u 2 
Au 3 (u - l) 2 



2GMn 



+21n(^ )+0(r- 2 ^ 



For comparison, we institute a sharper cutoff, 



Pt(x) 



M n 



T 



47rr 3 x(l + x) 2 (r 2 + x 2 ) 2 
The total mass of this profile is 



(A.25) 



(A.26) 



M Mr, 



T 



2(r 2 + If 



[2t 2 (t 2 - 3) lnr - (3r 2 - l)(r 2 + 1 - rvr)] . (A.27) 
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The projected mass density is 



r 2 47r(r 2 + l) E 



2(r 2 4- 1) 

[1 _ F{x)] + 8 F(x) 



+ 



+ 



r 4 -l 7r[4(r 2 + x 2 ) + r 2 + l] 



r 2 (r 2 + x 2 ) (r 2 + x 2 ) 3 /2 

r 2 (r 4 - 1) + (r 2 + x 2 )(3r 4 - 6r 2 - 1) 



r 3 (r 2 + x 2 ) 3 / 2 
The total projected mass is obtained as before, 



L(x) . 



M proj (x) = M 2(r2 ^ 1)3 [2[r 2 + 1 + A(x 2 - l)]F(x) 
+ - [vr(3r 2 - 1) + 2r(r 2 - 3) lnr] 



+ 



X/ 



[_ r 3 7r[4(r 2 +;r 2 ) _ r 2_ 1] 



+ [-r 2 (r 4 - 1) + (r 2 + x 2 )(3r 4 - 6r 2 - 1)] L(x)] }. 



Finally, the lensing potential is 

2GM 



-Wr 2 + it + 2(3r 4 - 6r 2 - l)rVr 2 + u L{y/u) 
+ 2r 4 (r 2 - 3) L 2 (v^) + 16r 4 ( M - l)F(v^) 

+ 2tV -3)(cos--L) 2 

+ r 2 [2r 2 (r 2 - 3) lnr - 3r 4 - 2r 2 + l] lnw 
+ 2r 2 r 2 (4rvr + (r 2 - 3) In 2 2 + 8 ln 2) 

- ln(2r)(l + 6r 2 - 3r 4 + r 2 (r 2 - 3) ln(2r)) 

- rvr(3r 2 - 1) ln(2r) 



As before, 



2G M proj (y^) 

,2 



i>\u) 

u 

The second derivative follows, 



2GM 



4(r : 



J ^y&{ [-8(«-l) + 2(1 -3r 2 )] F(V«) 
2(r 2 + 1) 



+ 



4(r 2 - 3) lnr + v ' \l - F(y/u)] + 8 
u — 1 

(r 2 + M )(3r 4 -6r 2 -l)-r 2 (r 4 -l) 
r 2 (r 2 + «) 



2w(3r 4 - 6r 2 - 1) r . „ 

+ VtT u/2 L v 7 ^ 
r 3 (r 2 + -u) 1 / 2 
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(3u + 2r 2 )[(r 2 + « 2 )(3r 4 - 6r 2 - 1) - r 2 (r 4 - 1)] 



r 3 (r 2 + w) 3 / 2 



L(y/U) 



- -(3r 2 - 1) 



r 



+ 



7T 



(r 2 + m) 3 /2 
Finally, the third derivative, 



2GM n 



T 



[2(r 2 + m)(3t 2 - 1) + u(4(r 2 + it) - r 2 - 1)}^.32) 
6(r 2 + 1) 20r 2 + 12 



c 2 8(r 2 + 1)% 3 
- 18 + 24(u-l) 



(u-iy 



u- 1 



30r 2 



fr 2 + w) 2 



6(r 2 + 1) ; 
(u- l) 2 



+ 



2(r 2 + l) 2 (5r 2 + 7) , 2 x2 , 2 x 
— + 4(r 2 + l) 2 (w + r 2 + 2) 



it — 1 

2\„.2 



+ (2 + 3r> 

,2 



(r 2 + 1) 
u- 1 



32 + 16(r 2 -3)lnr 



+ 



u 



t 2 (t 2 + u) 



+ -(3r 2 -l) 



4 [(r 2 + u)(3t 4 - 6r 2 - 1) - r 2 (r 4 - 1)] 
t 2 (t 2 + w) 

vr [24r 6 + 3m 2 (4u - 5) + 5r 2 w(9w - 4) + r 4 (60w - 8)] 
(r 2 + w) 5 / 2 

16r 10 + 30t 6 m(m - 4) + 9r 4 w 2 (w - 10) r . „ 
H ^r- — — Mv") 



+ 



r 3( r 2 + u )5/2 

-3(6r 2 + l)w 3 + 8t 8 (5m - 6) 



(A.33) 



r 3( r 2 +M )5/2 

All of these results reduce to the pure NFW case in the limit r — ► oo. 

As a final comparison, we reproduce the 0] formulas for a sharp cutoff at x = r. 
In addition, we will derive the third derivative of the potential. First, we define the 
auxiliary function 

1 



T{x) 



_j \Jt 2 — x 2 _ 1 \/t 2 — X 2 

tan — — — tan 



Vx 2 ^! 



ry/x 2 — 1 



(A.34) 



choosing the branch where T(x) is positive. This function is well defined everywhere 
provided we take the limit at x — 1, finding 



T(l) 



r- 1 
r + 1 



(A.35) 



Here, ArcTan from Mathematica and catan from C99 agree on the branch of tan 1 (x) 
for purely imaginary argument. 
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M proj (x) 



ln(l + r) 



T 



1 + r 



— x z 



r+1 



tanh 



Note that the total mass is 



M = M proj (r) 



M 



r z — ar 



ln(l + r) 



T(x) 



1 + r 



(A.37) 



(A.38) 



Note that these are valid only for x < r. For i > r, we simply have = and 

M proj (x) = M. 

We find that the potential involves at least the polylogarithm Li 2 ; we do not find a 
simple expression. We start with the first derivative, 

2G M proj (v^) 



Assuming u < r 2 , 



a 



(A.39) 



2GM 1 



c 2 it 2 



ln(l + r) 



r 



+ 



T — It 



1+r r+1 



— tanh 
2GM 2 

C 2 M 3 



w\/r 2 — u 3u — 2 



+ 



2(r +!)(«-!) 2(«-l 



T(v^)kA.40) 



ln(l + r) - 



1 + r 



+ 



r — u 
r+1 



— tanh 



.j \Jt 2 — u 



u^/t 2 — II 



2(r +!)(«- 1) 



15u 2 -20m + 8 , „ 

+ — ~ -^T(v^) 



it 



+ u[(u - l)r + u{2 + u)- r 2 (l + 2m)] j 



(A.41) 



8(l + r)Vr 2 l) 2 
Note that as u — > r 2 , diverges (see also Figure [3]). This is clearly because the surface 
mass density of this model is not smooth at u = r 2 . On the other hand, for u > r 2 , 
the derivatives can easily be computed as ip'(u) = 2GM/c 2 u, ip"(u) = —2GM/c 2 u 2 , and 
i)"\u) = 4GM/c 2 u 3 . 



Appendix B. The Sersic Profile 

The Sersic profile describes the light distribution of elliptical galaxies, and it also can 
be made to fit the mass distribution of haloes (33| . This profile is defined in projection, 

M 



rQir(2n) 



■ exp (— x 1//n ) 



(B.l) 



where r is a scale radius and x = r/r . The effective radius r e which contains half of 
the projected mass (or light) can be determined numerically. For the de Vaucouleurs 
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profile with n = 4, r e = 3459. 5r . The three-dimensional density distribution can be 
obtained by Abel inversion, 



p(x) = — 



dx 1 dY, 



vrr J x y/x' 2 - x 2 dx' 



tt/2 



/ du sec m - — (x' = x sec u 

7rr Jo 



(B.2) 



and we find in practice that the integral with finite range works well numerically As 
before, we need the total projected mass inside radius x, 

T(2n,x^ n ) 



M proj (x) = M I 1- 



T(2n) 



The total mass is 

M = M . 



(B.3) 



(B-4) 



The lensing potential can be expressed as a generalized hypergeometric function 351 ] . 

^(u) = 2G f°^„ 2 F 2 (2n, 2n; 2n + 1, 2n + 1; -u 1/2n ) , (B.5) 



c 2 (2n)! z zv ' ' ' ' 

where again u = x 2 . When 2n takes integer values, these results can be expressed in 
simpler terms as follows, 



if \ 2GM °o 
ip(u) = 5 — 2n 



+ exp(-n 1 / 2 ") 



l nu V2n + 7 _ Ei (_ n l/2n) 
j/2n- 



2n-2 



U J 



3=0 



with 



2n-l 



k=j+i 



(B.6) 



(B.7) 



and Ei(x) being the exponential integral function. Now we calculate the derivatives of 
ip. As before, 

2G M proj (v^) 



When 2n is an integer, we can express this in terms of elementary functions, 



2GM 1 



c 2 u 



j=0 



The second derivative is 

2GM 



u- 



l_ exp (-n 1 / 2 ™) J2 



j/2n ' 



(B.8) 



(B.9) 



i=o 



+ 



u(2n)! 



exp (-n 1 / 2 ") 



(B.10) 
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Finally, the third derivative is 

2GM 



c 2 



2 / , ito'^u^' 



V i=o J ' 

1 exp (-„'/-) ( 2 + ^ 



■u 2 (2n 

Appendix C. Elliptical Potentials 



(B.11) 



The power of considering ^ as a function of u = x 2 and not x becomes apparent when 
considering elliptical potentials. We will describe how the various observable quantities 
are obtained from the potential, and how these are modified for the case of an elliptical 
potential. 

The deflection angle is just the gradient of the potential, 

a = V*p(u) = ip'(u)Vu. (C.l) 
The lens equation is given by 

s = f-a(f), (C.2) 
and thus the magnification matrix is 

J = I- V®W>«). (C.3) 
For this we need the second derivatives of ip, 

V <g> Vip(u) = ip'(u)V <g> Vw + f (m)V« <g> Vm. (C.4) 

Now we can define u so that the isopotential lines are ellipses and not circles: 

u = (1 - e)x 2 + (1 + e)y 2 , (C.5) 

v = -(l-e)x 2 + (l + e)y 2 , (C.6) 

* -(£;£)■ 

and higher derivatives, such as the third needed for flexion, vanish. Deviations from 
ellipticity manifest as dependence on the variable v. 

It is well known that elliptical isopotentials can lead to unphysical projected mass 
densities. For the singular isothermal sphere, e > 0.2 gives isodensity contours that 
become peanut-shaped. We propose a simple resolution to this problem that will be 
adequate for most purposes: add several elliptical potentials at the same location. We 
find that three potentials, one with ellipticity e, one with ellipticity fe, and one circular 
potential, can be summed in such a way as to give nearly elliptical isodensity contours 
with ellipticities as large as e = 0.8, in other words axis ratios as large as three to one. 
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Consider a power-law density profile p{r) oc r 7 , yielding a potential ipiu) oc u^ 3 . 
We combine three such potentials, 

ip(u, v) = a e ip e (u) + a fe ip fe (u) + a ^o(w), (C.9) 

where the subscript indicates the ellipticity. Note that all three have purely elliptical 
isopotentials, but ellipses of one ellipticity are expressed with both the u and v variables 
of a different ellipticity Expanding in e, the v dependence can be canceled up through 
order e 2 with the following choice of coefficients, 

af - = /[(5-7) 2 -/P- 7) r (c ' 10) 

ao =(l-/)(|^/-lW (C.ll) 
A fraction / = 1/2 produces the best results in most cases, 

a e /2 = ^a e , (C.12) 
1 + 7 2(1 + 7) 

For potentials more complex than a pure power law, the cancellation is only exact at 
a single radius and hence the construction of the elliptical model which fits the wide 
range in radii is not trivial. We leave this issue for future work. 

It is also well known that elliptical isopotentials can lead to an inferred surface 
mass density that is negative. This occurs for any finite ellipticity when the asymptotic 
density profile is r -3 or steeper. The total inferred mass is always positive however, 
and the negative surface mass densities inferred are small enough to not be a serious 
concern. 

We make a note of the total mass of elliptical profiles here. Assume a profile 
described by a lensing potential ip(u) with u defined as above. In general, the surface 
density (inferred from the convergence) is given by 

S oc^j' + ^j"{u + ev). (C.14) 

If this inferred surface density contains a finite total mass, the total mass of the lens as 
a function of the ellipticity of the isopotentials is given by 

M tot (e) = - 7 =L= M tot (e = 0). (C.15) 
More practically, we relate the following projected elliptical power-law mass density 

£ = ^ [(1 - e iso )x 2 + (1 + e iso )y*} (1 " 7)/2 , (C.16) 
with that derived by combining three elliptical isopotentials 

tp = b e [a e i) e + a e/2 ipe/2 + a ipo] , (C17) 

where 

= [ (1 - e ^ + (1 + ^|,-^ 

3-7 



REFERENCES 24 

We find that the following fitting forms connect these two models in the range 1.2 < 
7<2.9: 

e = e iso + [0-05(2.1 - 7 ) 2 + 0.257]4 4(L8 - 7)2+2 - 9 , (C.19) 
b e = 1 + 0.193(-0.9 + 7) L38 4o + 0.0121(0.1 + 7) 4 '° 6 4o- (C20) 

This relation breaks down rapidly for ej SO > 0.8, corresponding to an axis ratio of 3:1. 
Two isodensity contours match with < 0.6% level at ej SO < 0.6 and < 3% level at 



